Dynamics of linear polymers in random media 
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We study phenomenological scaling theories of the polymer dynamics in random media, employing 
the existing scaling theories of polymer chains and the percolation statistics. We investigate both 
the Rouse and the Zimm model for Brownian dynamics and estimate the diffusion constant of 
the center-of-mass of the chain in such disordered media. For internal dynamics of the chain, we 
estimate the dynamic exponents. We propose similar scaling theory for the reptation dynamics of 
the chain in the framework of Flory theory for the disordered medium. The modifications in the 
case of correlated disordered are also discussed. 



The conformational statistics of linear polymers and their 
dynamical properties in porous or disordered media are 
being investigated extensively far a long time [1-14]. The 
universal scaling properties of a linear polymer chain in 
a good solvent, which is ideally modelled in terms of the 
self- avoiding walk (SAW) on a regular lattice, are now 
well established fljjl^ . However the situation is far from 
clear in the domain of a polymer moving in a disordered 
medium. Knowledge about the effects of structural dis- 
order of the medium (or the lattice) on the polymer (or 
the SAW statistics) and on its dynamical properties are 
important both for understanding of the general influ- 
ence of disorder on the critical behaviour and also for 
applications. 

The effect of structural disorder of the lattice on the 
SAW statistics, in particular when disorder itself is not 
critical, is still under investigation. The initial indica- 
tion 0] of the instability of the pure SAW fixed point 
in the effective n-vector model, was latter shown to be 
inappropriate D in the SAW [n 0) hmit. However, 
it seems to be established by now |^-|^] that the proper 
considerations of the quenched structural disorder of the 
underlying lattice (fixing one or both ends of the SAW 
on the dilute lattice indeed induces a modified crit- 
ical behaviour for the static conformational statistics of 
SAWs on disordered lattices. The polymer size exponent 
ly is defined through the relation Ra ^ N'^ connecting the 
radius of gyration (end-to-end distance of the SAW) Rq 
of the polymer chain to the molecular weight (SAW chain 
length ) N. The disordered media are usually modelled 
as randomly diluted lattices with lattice site or bond con- 
centration p. It is found that i'd{= v{p < 1)) is greater 
than the pure SAW size exponent iy{= v{p = 1)). Of 
course, at the percolation threshold p = Pc, the under- 
lying lattice becomes a fractal with effective dimension- 
alities less than the lattice dimension d. The SAW size 
exponent VpX= ^{p = Pc)) is therefore clearly different 



|,y and in fact is larger than and i^d- The multi- fractal 
efifects of the percolation cluster on the SAW conforma- 
tional statistics at p = pc has also been investigated re- 
cently . As mentioned already, all these investigations 
attempt to describe the static properties of polymers on 
random lattices. 

Compared to the above investigations on the static 



properties of the polymer configurations on disordered 
lattices, not much is known about the effects of quenched 
lattice disorder on the dynamics of such linear polymers 
Q. Recently, of course, the scaling theory and simu- 
lations for the reptation dynamics of polymer chains in 
weakly disordered lattices (p ~ 1) has been investigated 
|Q,|o|. The dynamics of un-binding of polymers in a 
random media has also been studied and the con- 
strained polymer dynamics has also been investigated us- 
ing a variational method [ p^ . In this report, we propose 
scaling theories of polymer dynamics in such porous or 
disordered medium using the scaling theory for polymers 
Jl5| , p^ and the percolation statistics [TtI] . 

We first consider the Rouse model [[l5[ where the inter- 
actions are local. The Brownian dynamics of the centre 
of mass Rg can be written as ||l^ 
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where f{t) denotes the temperature (T) dependent the 
random stochastic noise. One thus gets 
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for the average fluctuation in Rq. This gives the dif- 
fusion constant D = ARq/I — ksT/N. In the node- 
link- blob model [0 describing the percolation cluster, 
the occupied sites or bonds form a super-lattice with lat- 
tice constant £,p ^ {p ~ Pc) where £,p is the percola- 
tion correlation length. In this model, however another 
length-scale, namely "chemical length" (along the lattice 
edge) scaling as Ip ^ {p — Pc)~''^ , appears due to the tor- 
tocity of the percolating paths. Moreover, there exists 
random blobs (multiply connected regions) and dangling 
ends. Assuming the possibility of Brownian diffusion on 
such super-lattices, we find that the diffusion constant D 
of a Gaussian chain can be written as (cf. ||l5|,|l6|). 
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where 
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Here the radius of gyration Rq is scaled appropriately 
on the percolation cluster (super-lattice aX p > Pc) by 
5p :Rg — > Rc/ip- It may be noted here that the chain- 
length N appears here in the above diffusion equation 
as mass and thus remains unchanged and should not be 
scaled with the chemical length Ip (see however the later 
discussion on the reptation dynamics). Clearly, the ex- 
ponent an is always positive and therefore D decreases 
to zero here as p ^ Pc- 

If one considers the long-ranged (hydrodynamic) in- 
teractions on top of the Rouse-like interaction, as in the 
Zimm- model, one gets (cf. p^) 
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using similar scaling of the appropriate variables. Here 
g{r) ^ [N/Rg'^] g{r/RG) denotes the pair correlation 
function with Rq scaled by and /i(r) ^ r^^ denotes 
the mobility function. We therefore get 
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Here also the centre of mass diffusion constant D de- 
creases to zero following the above power law, as p — s- pc- 
For the internal dynamics or local excitations of the 
polymer chain, one can estimate the dynamic exponent 
z following the same scaling argument [|l6| for the pure 
case. Writing the dynamical structure factor as 



D 



~g{ty-DRG'-Hl-l 
~g{{p-p,ykH). 



(7) 



Since the local excitations (for kRc ^ 1) are inde- 
pendent of TV, assuming D ^ {p — pc)°'N^I^{-w\ih a = 
2i'p,(3 = 1 for the Rouse model (Eq. (3)) and a — i^p, 
(3 = Vd (Eq. (6)) for the Zimm model, one gets 

z = 2+—, and = a + vJz - 2). (8) 

Hence z = 3 for all dimensions in the Zimm model, while 
z changes to 2 + 1 /vd for the disordered case in the Rouse 
model and its magnitude becomes equal to 4 for > 6. 



We have assumed in the above discussion, the possi- 
bilities of the Rouse or Zimm dynamics of the polymer 
chain in good solvents within the pores of the random 
media. This obviously necessitates the correlated per- 
colation picture of the porous media, and hence the ap- 
propriate modifications [|2| of the polymer statistics in 
such medium (see below). However, very near the perco- 
lation threshold, the polymer dynamics can perhaps be 
only reptation type along the constrained pore tubes. 
Following the earlier investigations the diffusion 

constant D for the polymer, making reptation dynamics 
in the disordered media, can be written as 
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where Rg ~ N"'', t ^ and AF{Rg) denotes the fluc- 
tuation in the polymer free energy due to the structural 
disorder of the medium. A simple Flory type estimate 
for the free energy F can be written, following Smaller 
a/ 1711, as 
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where i'{= 3/((i -|- 2)) is the SAW size exponent on the 
pure lattice and the second term, arising due to the 
random fluctuation of the density-density (excluded vol- 
ume) interaction, gives the conflgurational fluctuation 
AF. The most important point here is the rescaling of 
the polymer mass with the chemical length Ip. Until the 
above equation, we considered scaling of only the fluctu- 
ation ARg by the lattice constant £,p of the percolating 
super-lattice, and not of the polymer mass TV. In the 
above expression for the Flory free energy, on the other 
hand, one has to scale as well the chain length or mass N 
hy Ip Q . This is necessary here to cast the free energy in 
the disorderd media in the pure lattice form; or in other 
words, as in the free energy, both temperature and mass 
renormalisatins are allowed. The minimisation, with re- 
spect to R of the above free-energy gives the polymer size 



as 
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The above result may also be obtained by scaling the re- 
lation {Rg ~ N'^'') connecting the radius of gyration to 
the polymer chain length: Rg — > Rg/^p ^^nd N N/lp 
ITsf . The value of the exponent 6 is positive [0 for 
d ^ 2 and 3. This indicates that the polymer swells (for 
fixed N) as p ^ Pc- The possible negative value of 6 for 
d > A may indicate localisation of the polymer in higher- 
dimnensions (Gaussian chains) in the random media M. 
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The fluctuation AF can therefore be written as 

(12) 

It may be noted from the above expression that Vd > v 
for aU d < 4 and we assume Vd = v = 1/2 for d > 4. 
The values of the exponents (f) and x ^-re positive for 
d < 4. The above Eqs. (10) and (12) together describes 
the scaling behaviour of the diffusion constant {D) of the 
reptation motion for the polymer chain in the disordered 
medium. We find that for d < A, D decreases exponen- 
tially not only as the chain length TV increases but also 
as p — > Pc- At d = 4, X vanishes while <f> becomes slightly 
negative so that the diffusion constant D in (12) becomes 
independent of N and decreses with p approaching pc- 
This seems to be related to the possible localisation of 
polymers in higher-dimensions as discussed already. 

As mentioned already, these dynamical properties of 
polymers in porous media demand the consideration of 
the correlated randomness of the percolating cluster. 
Very recently, Blavasts'ka et al ||l^ have studied the 
static scaling properties of polymers on a d-dimensional 
random lattice with "correlated" randomness that has 
a power-law fall of the form l/r" for large separation r 
pof . Similar Flory-type estimate for the free energy F 
in the case of relevant long-range correlation {a < d) will 
be given by Eq. (10) where the spatial dimensionality d 
in the second term will be replaced by the parameter a 
denoting the range of interaction. For a > 4, this fluctu- 
ation term becomes irrelevant and thus we get the mean 
field critical behaviour. Using therefore the expression 
(11) for the Flory estimate of i^d for the non-trivial cases 
a < 4 (d ~ 4), we get 

2 1 ^ 

'^d = r— V - 7T + — > (13) 
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where S = 4 - a and ly = 1/2 + (4 - d)/12. Clearly, 
this size exponent Vd, obtained via Flory theory, does 
not match the renormalisation group result |l|l even m 
the first-order in 6. This is is similar to the failure of the 
Flory approximation for the exponent v in the first-order 
in (4— d) even in the case of pure lattice (with short-range 
interactions) [|l5|. It may however be noted that this re- 
sult for Vd is for the long-range case a < d, where the 
renormalisation group study in ||l^ failed to locate any 
stable fixed point. For a > d, the above method clearly 
yield the pure (short-range) SAW result. However, the 
renormalisation group study suggests the existence of a 
long-range fixed point in this region [ p^ . 

In conclusion, we have studied the dynamical proper- 
ties of lattice models of a linear polymer on a percolat- 
ing lattice model of the random or porous media, using 
simple scaling ideas for both the polymer chain and the 
disordered media. The results for the dynamical expo- 
nent z describing the internal dynamics, and the diffu- 
sion coefficient D of the center of mass of the polymer 



chain, using both the Rouse and the Zimm model, have 
been investigated for such disordered media. We find 
that while the diffusion coefficient decreases with p — pc, 
following power laws in both the cases, the dynamic ex- 
ponent z remains unchanged to its pure value 3 for the 
Zimm model but gets modified for the Rouse model. Ap- 
plications of the Flory theory for the polymer free-energy 
and the scaling theory of the percolation disorder indi- 
cate that the reptation diffusion coefficient D of the poly- 
mer chain decreases exponentially with the chain length 
(N) and also inversely with the deviation {p — pc) from 
the percolation threshold Pc- D ~ eyip[—N^{p — pc)'^], 
where the exponents x ^■iid are both positive and re- 
lated to the polymer size and the percolation exponents 
respectively. Thus, based on a simple Flory-type theory, 
we have extended a few dynamical properties of linear 
polymers near the percolation threshold and have also 
discussed the case of correlated randomness. 
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